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I. INTRODUCTION 

The free energy and the static correlation functions of a gas of identical particles with equally spaced energy levels 
can be calculated exactly using a combination of the path integral method the method of symmetrized 

density matrices ^ and inversion techniques for generating functions |^,^ . The free energy, the internal energy, the 
specific heat, the moments of inertia the density and the pair correlation function have been worked out earlier 
for confined bosons. For the fermion case, the free energy, the ground state energy and the energy of the Fermi level 
were also studied before Using the same formalism, we consider in the present paper the density and the pair 
correlation function of a spin polarized fermion gas with equally spaced energy levels. Most of the calculations are 
analytical, and numerical work is only required for making the graphs. Nevertheless great attention had to be paid 
to the accurate calculability of the expressions in view of the sign problem originating from the statistics. 



■1^ 

Inspired by the recently observed Bose- Einstein condensation t-lOL much theoretical work has been done on boson 
models with equally spaced energy levels using other methods |ll|-|l7|. Analogous models for fermions, taking into 
account the confinement as well as the statistics, have been proposed and are studied in (l|Jl|l as a model for a 
quantum dot, and in ||2^,|2l| as a model for confined fermion alkali metal vapors. 

The model that we have been using before and that we also will investigate in this paper consists of N identical 
particles with mass m in a harmonic one-body confinement potential given by: 

and interacting with each other trough a two-body potential given by: 

N 

V^. = -^E(r.-,f. (1.2) 

k> ' This model has been studied for distinguishable particles |Q where it turns out to be equivalent to a set of oscillators. 
; I ' Three oscillators have a frequency O and are related to the degrees of freedom of the center of mass; the remaining 
5^ ] internal degrees of freedom oscillate with a frequency w — — N^. The model has led to some controversy 

when put in an occupation number version to deal with the statistics of the particles p^ . Avoiding the occupation 
number representation by a direct projection on the symmetric or antisymmetric representations of the permutation 
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group Sn two things become clear: first the center-of-mass coordinate factorizes out, indicating independence of the 
internal degrees of freedom, and second the propagator associated with the center of mass remains invariant under 
projection, indicating that the evolution of the center-of-mass degrees of freedom of distinguishable or indistinguishable 
oscillators are the same. This aspect of the model may be clarified by the introduction of the center-of-mass coordinate 
R = J2jLi i^to two-body potential: 

1=1 

The rewritten two-body potential makes it clear that the center of mass oscillates with a frequency il = {N'y + w"^) 
that is lower or higher than the frequency w of the internal degrees of freedom depending on the sign of 7 , which 
we earlier denoted by uj'^ or — depending on the case. For repulsion 7 > there is a stability constraint on the 
confinement potential: H. has to be large enough to keep all the repelling particles together. 

Because this condition on Q depends on the number N of particles, an approach with a fixed number of particles 
is mandatory from the very beginning, resulting in a constraint on the summation over the cycles in the cyclic 
decomposition of the permutations. We circumvented the complications of this constraint by first transforming to 
generating functions, and subsequently inverting the transforms to obtain the partition function, the density and 
the pair correlation function. The application of these mathematical techniques are among the main results of this 
paper. The fact that they allow to derive explicit expressions for the density and the pair correlation function of an 
interacting fermion system is also a new result. 

The paper is organized as follows. In section II, we collect the expressions obtained before for the fermion case, 
and in section III we present the calculation of the density and the pair correlation function in a general theoretical 
setting. We show how to evaluate these response functions for a system with minimal internal energy and a given 
number of particles. In the last section we discuss the results, and put the model in perspective. 



II. FERMION OSCILLATORS 



In this section the basic formulas which have been derived before for identical particles (bosons or fermions) are 
summarized and rewritten in such a way that they are more appropriate for dealing with fermions, in particular in 
view of the numerical treatment. 



A. The recurrence relation approach 

In our path-integral treatment Q , a recurrence relation was obtained for the partition function Z/ (N) corresponding 
to the degrees of freedom with frequency w in the relative-coordinate system. Introducing b = e~^'^ for brevity in 
the notations, we found that 

= ]^ E T3pv^ ^' M • (2-1) 

rn=0 \ / 

This recurrence relation applies for bosons — +1) and for fermions — —1). The subscript / refers to identical 
particles, which can be specified to be fermions (subscript F) or bosons (subscript B). A similar recurrence technique 
was used in |^ to obtain the contribution of the relative-coordinate system to the spatial Fourier transform of the 
density 



1 exp (-^ coth ^l(3hw) (iV - /) 



The center-of-mass contribution introduces the following factor: 

For the Fourier transform of the pair correlation function we found that 



(2.3) 
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where 



(1 - fej) (1 - u-jy 



(2.5) 



The center of mass does not contribute directly to the pair correlation function. It was shown in Q that the actual 
numerical implementation for fermions of (2.1) suffers from a sign problem. Therefore a generating function approach, 
followed by an inversion of the transform, turned out to be a more appropriate method of evaluation. 



B. The generating function approach 

The generating function Sj? (/?, u) for the partition function 

oo 

~f{P,u)=Y,I.f{P\N)u^ (2.6) 

was obtained before ||^. It should however be noted that this function is a calculation tool to obtain the partition 
function "Lp {I3\ N) where the number N of particles is given and not subject to fluctuations. A discussion of this 
point can be found in ||24[| . This discussion points out what the difference is between an ensemble-based approach to 
the problem (see e.g. p5| ) and quantum statistics for a finite number of particles. The generating function for the 
Fourier transform of the density is: 

oo 

e„(ii,q)= ^Z;^(/3|7V)Ar^iq«A'. (2.7) 
The generating function for the Fourier transform of the pair correlation function is: 

oo 

gg(«,q)= ^ZF(/3|iV)iV(7V-l)gqU^. (2.8) 
The defining equations for the density and the pair correlation function 

"W =^(|:''('----')) =/^",«-''. (") 

are taken over from |^ , and will be rewritten in a numerically tractable form. In comparison to Q , the normalization 
factor in the defining equation for the pair correlation function has been changed from N to N {N — 1) , such that 
J cPrg (r) = 1. Considering a model of N fermions with parallel spin in an harmonic confinement potential, and with 
an quadratic interparticle interaction as discussed in the introduction, the following generating functions are obtained: 



^f{u) ^ (1-6^)' V 1 
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with 



= (2.14) 

This means that for a system of non-interacting oscillators with eigenfrequency w, 'E.p (u) is formally a "grand 
canonical partition" function. However, strictly speaking it is not the grand-canonical partition function of the 
system with interaction for two reasons: first one has to take the center-of-mass correction into account, and second 
the eigenfrequency w in the relative coordinate system depends on the number of particles. But given w the full 
mechanism of generating functions in the relative coordinate system is applicable, provided afterwards the necessary 
center-of-mass corrections are taken into account. 

C. The density and the pair correlation function 

From the Fourier transform of Qn (w, q) and Qg (u, q) , the density in the relative coordinate system and the pair 
correlation function can be obtained in real space: 



where 



h 

Including the center-of-mass correction, the Fourier transform of (|]^) has to be used instead of riq, giving 



(2.17) 



with 



Ai = , , . (2.19) 

coth ^/3?iu; 



D. The ground state correlations 

The expressions for the density and the pair correlation will be studied in the zero-temperature limit 6^0, i.e. in 
the ground state. The relation between the number of particles and the fugacity implies: 

^^2> 'l + uV>l'V ^ (1-b'f 

in which case the chemical potential ji in the fugacity u = e^^ becomes the Fermi energy. Integrating the generating 
functions as follows: 



r g» (<q) ^^/ =f2^F (N) nqu^ = (u) (nq) , (2.21) 

^0 " N=0 
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\ \ ^^f;^du''du'=J2zAN)9^u'^ = E^{u){g^). (2.22) 
Jo Jo iu") 



/o Jo {u"Y j^^Q 
the Fourier transform of the density and the pair correlation function become: 



' riHw.oftia^»,>f-.'i±i;)..', (2.3, 



~ (z.) Jo u'^^ -'^ (1 _ b^f V l-b' 



^f{u)Jo Jo (Wf hfy (1-60 

X ^exp (^-^) - (<3<!j)^exp {-2K^Qe^j)^ du"du'. (2.24) 

A straightforward but tedious calculation shows the mathematical equivalence between the expressions obtained in 
the recursion approach and those obtained by integrating the generating functions. For the density this calculation 
proceeds as follows: 
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The equivalence for the pair correlation function can be obtained along the same lines. 

III. THE DENSITY IN THE GROUND STATE 

In this section the ground state density will be calculated, thereby introducing the Fermi level explicitly. Simulta- 
neously, the sum rules will be checked. 

A. Fourier transform of the density at low temperature 



Starting from the generating function ( 2.12 ), the series expansion of the exponent in Qn (u, q) e!^ as a Taylor series 
in b gives 

g„ {u, q) e-' ^ ^ (-2^^)" ^ r (fc + 3 + n) 

Epiu) to ^r(fc + l)r(3 + n) i + 
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for which the analytic continuation has to be found as a function of the fugacity u = e '^'^ for 6 = e ^^"^ arbitrarily 
small. If L denotes the lowest unoccupied level, this expression becomes, using u&i = h-'^+°' (with a -> 0+): 

g„ {u, q) e-' _ ^ {~2n^Y ^ T (fc + 3 + n) 

n! ^ r(fc + l)r(3 + n) l + 

In the low-temperature limit (6 — > 0) it is clear that l)-L+a+n+k j (^^ + tends to zero for -L + n + k > 

and to unity if —L + n + k < 0. For T ^ the summations can then be restricted to n < i — 1 and k < L — 1 — n, 
and Qf.i {u, q) /Sp (u) becomes a polynomial in (— 2k^): 

g„(^,g)e-^ _^^ (-2^^)" ^^" r(fc + n + 3) ^ T (£ + 3) „ 

SfN ^„ n! A. r(fc + l)r(n + 3) T (L - n) T (n + 1) T (n + 4) ^ ^ 

We note in passing that this expression can also be written as ^L{L+l){L + 2) iFi (l — L;4;2k^) where iFi 
denotes the confluent hypergeometric function which is related to the generalized Laguerre polynomials _L„ " (z) — 
^^"''„"'"^ iFi {-n; l + a;z) where {a,n) = a{a + 1) ■ ■ ■ {a + n - 1) = HlZ^ (a + k) = ^^p°^"^ . 
For q = one readily finds Gn (u, 0) fE-F (u) = N. One then obtains from hq = Gi [u, q) /Gi (w, g = 0) in the limit 
0: 

~ - X (--!^\\^ 6r {L) /^„_^y 



The Fourier transform of the density including the center-of-mass correction readily follows from Eq. (2.2). It 
essentially modifies the spatial decay, not the polynomial which arises from the fermion statistics. Taking into 
account that coth ^fjMl — > 1 and coth ^phw 1 for /? ^ 00 it follows that 



with 



p f 6r(L) / hq'Y 

"""^^ 4mVFy ^ r(l + n)r(4 + n)r(L-n) ^ 2mw J ' 



B. The ground state density 

The probability density of the fermions can be obtained from the Fourier transform np (r) = j ^j^'^F.qe^*'' ''. 
Using ( ^.5[ ), the angular integrations are readily performed, and the remaining radial integrations are of the form 
/~ Ki+2"e-''' sin {2Ka) dn = ^e-'" Ct"") a^+^"-^'T {j + i) . The duplication rule {2j + 1) = 

(j + r (j + 1) allows to rewrite T (j + 5) in gamma functions with integer argument. After some algebraic 
manipulations one is left with 

np (r) f^) e--"^ V 6r(£)r(2 + 2j) / 2I^y" 

^^^"UW ^r(i + ,)r(4 + ,)r(L-j)l 4u;j " 

t ... .',.^.J -^^- (3.7) 



^^r(j-n+l)r(2 + 2n) 



w 



The remaining double sum represents a polynomial in p^, the numerical evaluation of which presents no difficulties. 
The probability density in the origin r = becomes 

„ mi /mW^y/^^Y XWV 6r(£)r(2 + 2j) 

= j ^„ r 2^ j n(i + .)r(4 + ,)r(L-,r ^'-'^ 

in which the remaining sum can be identified to be the hypergeometric series 2 Fi (|,l-i;4;2.^), although this 
knowledge still requires further numerical treatment. 
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1. Sum rule for the density 



One can check that the density satisfies the relation J drnp (r) = 1 as it should. For that purpose, performing the 
angular integrations one is left with 

, , ,_24^7i^y r(L)r(2 + 2j) ^ (-4)" 
^""^^"^-0F^^lv r(i + j)r(4 + j)r(L-,);^^ro- + i-n)r(2 + 2n)yo ' ^ ^ 

= 6 V Li^^V r(L)r(2 + 2i) ^ (-1)" 

2wJ r{i + j)r(4 + j)r{L-j):^^r{-n + j + i)r(i + n) 

For J = 0, the summation over n yields 1, whereas for j > the summation over n gives the binomial series 
(1 - 1)^ = 0. Therefore 



f ( \WV 



r(L)r(2 + 2j) 



+ j)r(4 + j)r(L-j) 



= 1 



and the required sum rule is indeed satisfied. 



2. The mean square distance 

A analogous calculation allows to calculate the mean square distance to the origin, {r^) p ~ J drr^np (r) , and to 
derive a scaling factor \/{r^)p appropriate for distances. Similarly as for the sum rule above one finds 

3/^ l^y r(£)r(2 + 2j) ' 3 + 2n 



mWf^^y 2wJ T{l+j)T{A + j)TiL-j)^^' F (j + 1 - n) F (1 + n) ' 
Separating out the contribution of j = to the summation, one is left with 

/^2\ 3ft (l l^y F(L)F(2 + 2j) ' 3 + 2n 

mW ^2 f^^y 2w) F(l+j)r(4 + j)F(L-j),^/ ' F (j + 1 - n) F (1 + 7i) 

Because I]^=o (~^)" r(n+i)r(-n+j+i) = r(i+j) ~ ~ r§) ^ ^f^^ ^' 1™^* ^ 1 can easily be taken: 

3ft 3ft, ^ 3ftu;-M^ 3ft 

\ (L - 1) = 1 

2 mW 4 mw 2 m Ww 4 mw 

Hence the mean square distance of the spin-polarized fermions from the origin is proportional to L N^/'^ for iV 3> 1. 



(r^)^^ \ L-1 = 1 [L + l). (3.9) 



3. Numerical results 

Having established the relevant distance scale, the scaled density np (r) /np (0) is plotted in Fig. 1 (for L = 1 till 
5) and Fig. 2 (for L = 6 till 10) as a function of p = r^/ mw/h for the case of non- interacting fermions {W = — w). 
The density in the origin for these cases is given in Table 1 . The density profiles can be compared with those obtained 
in ref. pT[ | because the canonical ensemble and the quantum statistical partition function lead to the same predictions 
for this model provided the number of particles is given. Actually, we did not make a detailed comparison because our 
calculation is done for the ground state while the calculations of are done for finite temperature. Nevertheless the 
results for completely filled Fermi levels look very similar. The plots can also be made for the interacting case Q w, 
but the differences are minor. For repulsive interactions the condition < w < has to be taken into account. From 



the defining equation (3.6) for W this leads to 



N W , , 

and therefore the center-of-mass contribution to the density is washed out by the other degrees of freedom, except for 
a very limited number of fermions. 
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IV. THE PAIR CORRELATION FUNCTION IN THE GROUND STATE 



In this section we will repeat the same analysis as in the preceding section, but for the pair correlation function. 

A. The ground state expressions 

In contrast to the calculation of the density, which was obtained from its Fourier transform, the calculation of the 
pair correlation function is easier in real space, starting from ( ^.16 ). Because of the decay proportional to e^^'' , we 



expect that gp (r) e^P is a polyno mial i n r, which we want to be able to compute with sufficient accuracy in the 
low-temperature limit 6^0. From ( ^.16| ) one readily obtains 

" TV (TV - 1) ^ ^ (1 _ 6^-)3/2 (1 _ yf^ 

X ^ ^ '-^ . (4.1 

(1 - b^+if/^ 

Expanding the exponential in a power series, and introducing the Taylor series for Qj+ij ( ^.Ej ) and for the denom- 
inators, the resulting 5-tuple series expansion can be rearranged in such a way that it can be summed in the limit 
& — > 0. The mathematical details of this derivation are given in appendix A. The following expression was obtained 
for the low-temperature limit: 

^ e-2P /mw\ 3/2 (p72) 

^^^'■^ NiN-i)\27rh) A.Z.Z.r(fc + i)r(? + i)r(L-j)r(i + j-or(i + j-fc) 



(4.2) 



The behavior in the origin is determined from A; = and I = 0, and obviously results in the expected result 
gp (0) = 0. In Fig. 3 the pair correlation function of the spin-polarized fermions is shown as a function of p for various 
Fermi levels. 

B. Sum rule for the pair correlation function 

An important check on the validity of the derivation above is provided by the condition that the pair cor- 
relation function has to satisfy the sum rule Jdrgplr) — 1. Performing the angular integrations, and using 

/~ p2 (p2/2)'^'+" e-^p"dp = V2r (I + n + k) one obtains 

2 T^{l+j-k-l)T{\+L-j + k + l) 

''^^^'■^ - Tv(iv-i)y^^^^^^^r(fc-n)r(?-i-i)r(L-j)r(i + j-or(-fc + i + j)r(|-fc)r(|-/) 
2 g^^. (-i)'+'=rM| + ,)r(| + / + fc)r(i + L-,) 



iv(A^-i)0F^^^^^^r(fc + i)r(/ + i)r(i-j)r(i-fj-or(-fc + i+j)r(|-KOr(| + A:)- 

We checked analytically for i = 2, 3, ... 10, that these summations indeed yield 1, and therefore are rather confident 
that it holds in general. 

C. The Fourier transform 

Although it is possible to calculate the Fourier transform oi gp (r) , this calculation is rather involved, and it is 
easier to evaluate the zero-temperature limit of directly from the generating function along the same lines as above. 
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The actual calculation proceeds along the same lines as for the density, the Fourier transformation of the density, 
and the pair correlation function itself: expand the expressions in powers of 6, and study the analytic continuation of 
Qg (it, q) in the fugacity u — e*^^ for ub^^^ — The calculation only differs in the details from the three examples 

in the previous subsections and it seems pointless to report the full details of the calculation. Again using ~ ' 
the result is 

1 ^ 1 {-2^T r(2~j+n + L) 



r (L - j) mV (-m + 1 + j) ^^^^^ r2 [n - j) {n ~ injT [i + n)T [j - n + m + 1) 

One readily checks that this expression indeed has the correct long wavelength limit 

^'^""'^^ - ^\l{L + 1){L + 2)\ -\l{L + 1){L + 2) = N{N~1). (4.3) 



{u) \6 / 6 



With the normalization implied by the defining equation ( 2.10| ), one finds the following polynomial in from 
c/q = Qg ("", q) /Gg (u, q = 0) 



2 



"U T{L) +;^^r(l + n)r(L-n)r(4 + n)^ '"M 6 T (L) 



vf H^^)"^^' r(2 + n + L) ^ I 

^ r2 (n) r (L - j) r (3 + n + j) {n + j - m) mT {j - m + 1)T {m + 1 - n)' ^ ' ' 

Further simplifications are not easy to obtain, but this expression does not present numerical problems. 

V. CONCLUSION AND DISCUSSION 

In this paper we have shown that the density and the pair correlation function of the spin polarized fermion oscillator 
model can be calculated taking exactly the statistics of the particles into account. We studied explicitly the ground 
state correlations of this confined fermion system for a set of particle numbers that correspond with fully occupied 
Fermi levels. One of the reason for this limitation is that there is an additional simplification in the algebra, which is 
already tedious. Another reason is that this choice allows one to compare results calculated by inverting generating 
functions with the results based on the iteration of the partition function. For a limited number of Fermi levels 
{L — 2, 3, 4) both methods have been used to check the results depicted in the figures. 

The wiggles of the density as a function of the distance from the center are more pronounced with an increasing 
number of particles and are therefore clearly related to the increasing density of states at the Fermi level. Also the 
trend that the most probable distance between a pair of fermions decreases while the range of distances wherever 
the pairs can be found increases with an increasing number of particles is worthwhile to remark. The mean distance 
and its variance are also calculated, because these quantities also provide important averages needed to perform a 
Jensen- Fey nman variational calculation for a confined system of fermions with a realistic interparticle potential and/or 
a more complex confining potential. The model has also some importance in itself, because it can be used to test new 
approaches to Monte Carlo simulations of interacting fermions such as many-body diffusion [p6 27 2^. 



In summary we have been able to calculate the thermodynamics and the static response functions of the harmonic 
model for a finite number of fermions without any recourse on the thermodynamical limit or on the theory of ensembles. 
Our approach is fully quantum mechanical and relies on the probability assignment for quantum systems in equilibrium 
and as required by the statistics on the projection on the symmetric (bosons) or antisymmetric representation of 
the permutation group. Of course the generating function technique borrowed from the mathematics of stochastic 
processes together with path integral methods are essential in the formulation as well as in the evaluation of the 
thermodynamical quantities and the static response functions. 
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APPENDIX A: MATHEMATICAL DETAILS FOR THE PAIR CORRELATION FUNCTION 



In this appen dix, a possible way to derive the low-temperature limit (12) from (|4.lD is given. We first expand the 
exponentials in ([4.1[) in a power series, followed by filling out Qj+ij. This leads to 



/ \ 1 -,oo_,/_, \noocxD ... 

<7F(r)e2P _ 1 1 /I \ ^,3/2\^+' 



. 27r rJ "=0 1=1 J = l 



1 (-ir 



^(1 - b^f^^ (1 - ¥f^^ (1 - (1 - 6')"+^/^ (1 - 6^)"+^/^ (1 - bi+'f^^ 



Using the binomial series [6' (l - b>) + b> (l - 5')]" = Y2=q Qb^^W'''^-'''^ (l - V)'' (l - 6')" ^ and the Taylor 
series expansion 

for 6 — > 0, one obtains the following quite involved series expansion 

/ r (§ -I- m - fc) F (I + p - n + fc) F (I + g + n) (~1)" F (| -|- n - fc -f m) F (| + fc + p) F (| + g) \ 
^ 1^ F(i-fc)F(f-n + fc)F(|+n) F (| + n - fc) F (| + fc) F (|) ) 

The summations over I and j can readily be performed, and using ub^/"^ ~ j^-L+a ^j-^^g result can be written as 
fff (r)e^p' _ 1 ^^^^^ (pV2)'+" 5-L+«+»+m+g 

( F (I - fc -H m) F (I - n + p) F (f + fc + n + g) _ (-1)^+" F (| + n + m) F (f + fc + p) F (f + g) \ 
V F(f-fc)F(|-n)F(f + fc + n) r (§ + n) F (§ + fc) F (|) J ' 

The analysis of the zero-temperature limit 6^0 then proceeds similarly as for the density. The fraction 
&--^+"+'=+P+V (1 -h &-i+"+fe+P+9) tends to unity if fc -h p -h g < L - 1 and to zero otherwise. This means that 
the summations can be restricted tofc<L— 1 — g, p<_L — 1 — g — fc and g < L — 1. A similar analysis applies to 
fo-^^+^+^+^+V (1 -I- 5-i+"+»+™+9)^ leaving a polynomial in pV2 for T ^ 0: 

gpjrje^P ^ 1 V V V V 

^'m^«^3/2 ^ N 2^ 2^ 2^ fc!?i!r7i!n!o! 

l2?ly 9=0 fc=0 n=0 m=0 p=0 ^ ^ 

/ F(|-fc + m)F(|-n + p)F(f + fc + n + g) _ F (| + n + m) F (| + fc + p) F (f + g) \ 

1^ r(i-fc)F(|-n)F(| + fc + n) ^ ^ T (| + n) F (| + fc) F (|) J' 

Also the summations over p an m can be done analytically. These are of the form 
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^ r (fc + I + p) _ 2fc + 3 + 2M 2-i-2fc-2M ^ (2fc + 2M + 2) _ F (fc + + |) 
^ p\ ~ 2fc + 3 M!(A: + M)! ~ (A;+|)M! ' 

and hence after some manipulations with the summation indices: 

^ e 2P /mwN 3/2 (pV2) 

^^^"■^ N \2TTh) ^^^^^r(fc + i)r(/ + i)r(L-i)r(i + j-or(i + j-fc) 



'r2(|+j-fc-0r(i + i-j + fc + (-ir'r2(|+j)r(i + L- 



r(|-fc)r(|-Or(| + A; + r(| + Or(| + A;)r(|) 
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L N (-rrh/mwf^^ np (0) 

11 1 

2 4 1/4 

3 10 1/4 

4 20 1/8 

5 35 1/8 

6 56 5/64 

7 84 5/64 

8 120 7/128 

9 165 7/128 

10 220 21/512 



TABLE I. Number of particles N and density in the ori- 
gin for the lowest unoccupied energy levels characterized by 
L = 1, . . . , 10. 



Figure captions 

Fig. 1: Scaled density np {'r) /np (0) as a function of 

p = ry^mw/h for the lowest unoccupied levels char- 
acterized by L = 1, 2, 3, 4, 5. 

Fig. 2: Scaled density np {'r) /np (0) as a function of 
p = r^mw/h for the lowest unoccupied levels char- 
acterized by -L = 6, 7, 8, 9, 10. 

Fig. 3: Reduced pair correlation function 7(r) = 

gp (r) {Trh/2mw)^^^ as a function of p = rsjmwjh, 
for the lowest unoccupied levels characterized by 
L = 2,3, ...,10, indicated by the numbers in the 
corresponding curves. 
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